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In a setup where two ferromagnetic electrodes are attached to a superconductor, Andreev bound 
states are induced at both ferromagnet/superconductor (FM/SC) interfaces. We study how these 
states propagate through the SC and interact with each other. As a result of this interaction, the 
energetic positions of the Andreev states are not anymore determined solely by the magnetic prop- 
erties of a single interface, but also depend on the interface distance and the relative magnetization 
orientation of the FM contacts. We discuss how these bound states show up as distinct peaks in the 
nonlocal conductance signal, and lead to marked asymmetries with respect to the applied voltage. 
We relate our results to nonlocal crossed Andreev and elastic co-tunneling processes. 
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The presence of Andreev bound states (ABS) at en- 
ergies below the superconducting gap plays a prominent 
role in transport through heterostructures involving su- 
perconductors, and forms a topic of continuing and long- 
standing interest. Such states can arise, e.g., due to 
multiple Andreev reflections at both interfaces of a SC- 
insulator-SC junction, and are directly related to the 
Josephson current through the junction [1] . Furthermore, 
the importance of a magnetically active interface between 
a SC and FM has been pointed out. Scattering on such 
a surface leads to spin-mixing and the creation of subgap 
ABS at the SC/FM interface 0, [sj. Such states influ- 
ence theproperties of, e.g., single SC/FM tunnel junc- 
tions |3,[5| and SC/FM/SC Josephson junctions 

Andreev reflection processes assume a nonlocal (NL) 
form when two metallic leads are connected to the SC: 
the incoming quasiparticle in one lead can be reflected as 
a quasihole in the other lead [8], giving rise to a nega- 
tive NL conductance. This so-called crossed Andreev re- 
flection (CAR) competes with elastic co-tunnehng (EC) 
(and, at high temperatures, charge imbalance effects [9]), 
in which the quasiparticle is transferred from one lead 
to the other via a virtual subgap state in the SC and 
gives a positive contribution to the NL conductance [lo| . 
CAR and EC processes cancel each other for tunnel con- 
tacts pjj. This cancelation is lifted at hi gher orders in 
the transmission [l2[, for FM contacts 11 -l^? oi" in the 
presence of interactions [l5|. Disorder effects have been 
addressed in Refs. [l6|. Even though solid theoretical 
progress has been made, the role of Andreev bound states 
in the nonlocal conductance has not been elucidated. 

In this Letter, we consider the non-local setup shown 
in Fig. [TJ where two FM point contacts are attached to a 
ballistic SC region. At each of the two SC/FM contacts 
subgap ABS form, which show up in the Andreev spec- 
trum [sl, [HI . They propagate on a coherence length scale ^ 
through the SC. Interestingly, these states interact with 



each other so that their energetic position is not only 
dependent on the magnetic properties of the separate in- 
terfaces, but also on the distance between the interfaces 
and their relative magnetization orientation. We study 
the profound influence of this ABS interaction on the NL 
conductance through the device when voltages and 
are applied across the contacts. For identical contact pa- 
rameters, we flnd that for parallel (P) magnetization, the 
NL conductance dl^i/dVi, is asymmetric in l^,, whereas 
it stays symmetric in the antiparallel (AP) conflguration. 
Such (a) symmetries will be explained with an intuitive 
picture based on the ABS positions in the system. 

Transport properties of heterostructures involving SCs 
are conveniently described using quasiclassical Green's 
functions (QC GFs), g{vY-,R^s)^ which depend on the 
Fermi velocity iTp, the spatial coordinate R and the quasi- 
particle energy e, and obey the Eilenberger equation 



\hvY ' V^g + [eh - A, ^] = 6, 
subject to the normalization condition g'^ 



(1) 
The 



" hat" refers to the 2x2 matrix structure of the propaga- 
tor in particle- hole (Nambu) space, and A is the SC order 
parameter. Within our approach, the exchange energy is 
incorporated by different Fermi velocities vyo- and mo- 
menta pfct in the spin bands 0, Thus, all elements 
of the SC propagators, ^sc, are 2x2 spin- matrices, while 
the elements of the FM propagators, ^cr, are scalars. 
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FIG. 1: (Color online) Sketch of the considered FM/SC/FM 
setup with two tunneling point contacts. A nonlocal scatter- 
ing event which involves both contacts is indicated. 
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FIG. 2: (Color online) Spin-up (positive axis) and spin-down 
(negative axis) DOS at different positions x for parallel (left 
column) and antiparallel magnetization (right column), (a) 
and (b): p = 0.9, L = ^, i^r = di^L = 0.47r. (c) and (d): 
p = 0.9, L = 2^, / ^R, with i^L = ±0.47r, = O.Ttt. 



A SC/FM interface enters the QC theory as bound- 
ary conditions relating outgoing and incoming propaga- 
tors . They are expressed in terms of the normal state 
scattering matrix S of the interface. We assume that the 
transmission (t"^) and reflection coefficients (r^) of the 
^-matrix are spin-diagonal with respect to the magne- 
tization of the FM, with diagonal elements exp(ii^^^) 
and r^exp(ii9^) where a G {t,i} and a, (3 e {SC,FM}. 

First we calculate the local density of states (DOS) at 
points along the line connecting both FM contacts and 
investigate the ABS interaction. We concentrate on the 
contribution arising from the trajectory that connects the 
two contacts (see Fig.[T]). Denoting x as the distance from 
the left electrode, we obtain from the solution of Eq. ^ 



-Re 



(i-rLafR^)/(i + rLafR^) , (2) 



where A^f is the normal state DOS in the SC, and with 
the notation a =t (|) when a =| (t). For subgap ener- 
gies, the Via- with i G {L, R} originating from the left (L) 
and right (R) contact, are given by [13. ITtI 

Tia{x,e) = [1^705^^ + i(7o£Si^ + A)tanh(l]xi/A^)] 

X [n-i{joAsi^^e)td.nh{nxi/A^)]~\ (3) 

and the Via result from Tia{x^ e) = r*^(x, —e). Here, Q = 

xi, = x^ xr = L — X (where L is the distance between the 
contacts), and ^ = hv^^ / A. Interface properties enter 
the DOS only via the product p = r^r^ and the spin- 



mixing angle 



1^1^ at left and right interface. 



In the following, we assume for simplicity p = Pl = PR, 
but allow for differing '^^l and '^^r. In particular, i^r = 
= when speaking of P (AP) magnetization. 



In Fig. [21 the DOS is shown at different positions x 
between the interfaces. ABS appear as distinct spin- 
resolved subgap peaks. In general, for a positive spin- 
mixing angle, a spin- up state appears at positive en- 
ergy. For negative spin- mixing angles, spin- up and spin- 
down peak positions are interchanged, as expected from 
Eq. (|3]). Fig. [2] reflects this behavior. For the P case 
'^^L = '^^R > 0, the peaks appear at identical energies at 
the left {x = 0) and right {x = L) interface. For AP 
magnetization with '^^l = — '^^r < [Fig. Efb)] the spin- 
up (down) state is created at a negative (positive) energy 
at the left interface, where 'i^L < 0, but at a positive (neg- 
ative) energy for the right interface, where i9r > 0. The 
DOS peak height decreases upon propagation through 
the SC. For AP magnetization, this leads to two spin- up 
(and down) states at each interface: the smaller one is 
created at the opposite interface and is attenuated upon 
propagation through the SC. This attenuation is less pro- 
nounced in the P case, where it is countered by the prop- 
agating state from the other interface. If I'^^lI 7^ |'^r|, 
ABS have different energies at left and right interface, 
and their spin-character depends on the sign of the spin- 
mixing angles [compare Figs. [2fc) and (d)]. 

These states correspond to poles in Eq. ([2]). From 
equating 1 + rLcrPRa = we obtain in the tunneling 
limit (pl = PR = 1) 

cos(V^ + ^) cos(V; + ^) = e-^ sin(^) sin(^) 

(4) 

for the spin-up ABS energies at the right interface. Spin- 
down states are found from Nc^{x^e) = N^{x^ —e). The 
energy dependence in Eq. (j4]) is contained in the variable 
il) = arcsin(6/A). ABS positions thus not only depend 
on the spin-mixing angle, but also on the interface dis- 
tance L and the magnetization configuration of the FMs. 
For lower p, i.e., higher interface transmissions, the DOS 
peaks are broadened into resonances. 

In Fig. [3l the spin-up DOS at the right interface 
{x = L) is plotted as a function of the interface distance 
L and the spin-mixing angle "d. For the P configuration, 
one obtains a single spin- up ABS at = A cos(|) in the 
limit L/^ oo. This state splits in two as soon as the 
interface distance L is decreased. This can be interpreted 
as follows: the ABSs originating at both interfaces [see 
Fig.[2ja)], having the same energy in the limit L/^ oc, 
repel each other as the interfaces are brought closer. 
To first order, the energy correction to is given by 



fe^ = Asin"(f)e-T«^"^2i, assuming Se <C \A - e^\. So 
the repulsion strength increases with i9. When the inter- 
faces are close enough, the higher energy state is pushed 
out of the SC gap; close to the gap, its energy can be ap- 
proximated by A - 2 A [(L/0 - cot ( f )] V [1 + (L/0 ^] ^ The 
critical length at which the state enters the continuum is 
thus given by L = ^cot(^). At these small lengths, the 
lower-energy state goes like £^ + 2A sin^(|)(l - e~T ^^"^) 



3 




FIG. 3: (Color online) Spin-up DOS at the right interface vs 
L (upper panel) and il) (lower panel), both for P (a,c) and AP 
(b,d) magnetization. Dashed lines are analytical approxima- 
tions (see text). Inset: propagation of the split ABS between 
the interfaces. Parameters in Figs (a,b): p — 0.9, 'd — O.Ttt. 
In Figs. (c,d): p = 0.9, L = 2.0^. Inset: ^ = O.Ttt, L = 2.0^. 



with £^ = Acos('i^) its energy for L = 0. When > ^, a 
zero bias state exists for L = ^In | tan(|)|. The propaga- 
tion of ABS split peaks through the SC is sketched in the 
inset of Fig. [3fc); the weight of the higher energy state 
is smaller than that of the lower one and decays faster. 

In the AP case, there are two ABS at energies = 
±Acos(^) in the limit L/^ oo. In terms of Fig. [2jb), 
the negative energy state corresponds to the ABS that 
has propagated from the other interface. Its weight is 
therefore reduced (it is zero for L/^ oo). As the in- 
terfaces are brought closer to each other, this state gains 
in weight and both ABS repel each other. To first or- 
der, their energies are ±£^ ± |A tan |e"~ ^''''^ 2). They 
both enter the continuum at L = 0. 

The influence of the spin-mixing angle is depicted in 
Figs. [3jc) and (d). Increasing i9 moves all states deeper 
inside the gap. For the P case, the second ABS only 
appears for values i9 > 2arctan(L/^). In the AP case, the 
negative-energy state is weakened as it must propagate 
through the SC to reach the right interface. 

These Andreev states show up in the transport prop- 
erties of FM/SC/FM heterostructures, in particular in 
the NL conductance Gnl- Assuming a constant current 
density j over the contact surface area, the NL current 
/r in the right contact is: 



SlSr cos(6'l) cos(6'r) 
L2 



(5) 



where <Sl(r) is the area of the left (right) contact, L is 
again the distance between the contacts, and ^l(r) is the 
impact angle at the left (right) interface (see Fig.[T]). Gnl 
is given by Gnl = OIr / dV^ , where Vl is the left lead volt- 
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FIG. 4: (Color online) CAR and EC contributions for P (a-c) 
and AP (d-f) magnetization of the FM contacts. Transmission 
coefficients: (a-c): = t^, (d-f): = tj, with |t^|2 = 0.14, 



0.07. For other parameters see legends of (b) and (e) 



age. The current density j at a point R in the right elec- 
trode isj=e^^N^J |^Tr(r3i;Fa^f (^Fa,^,^:)), where 
No- is the normal state spin-up (down) DOS in the right 
FM. Only the trajectory connecting both contacts con- 
tributes to Gnl, and the relevant Keldysh GF can 
be calculated from a generalization of Eq. ([T]) and the 
boundary conditions at the interfaces [l7| . 

The clearest ABS signature occurs for high spin- mixing 
angles (ABS deep inside the gap), and for the tunneling 
limit (ABS narrow and well-defined). Fig. |4] shows Gnl 
for such a set of parameters. A peaked structure is visi- 
ble, hinting at the subgap states, with a marked voltage 
asymmetry in the case of P magnetization. These obser- 
vations will now be explained with the help of Fig. O 

For voltages eV^ that align the Fermi energy in the 
left FM with a spin-resolved ABS energy in the SG, 
a quasiparticle with the corresponding spin can tunnel 
from the FM into the SG. In case of EG, the particle 
has to tunnel to the same spin band at the same en- 
ergy in the right electrode, while for GAR, an electron 
with energy —e^ in the opposite spin band at the right 
interface is absorbed by the SG to form a Cooper pair 
(see Fig. [5j). Both GAR and EG thus involve (to lowest 
order) two transmission processes, and the magnitude of 
the corresponding conductance peak is given by -7^^ 
for EG and T^^ • Tj^ for GAR, where the T are effec- 
tive transmission coefficients involving the product of the 
relevant spin-resolved transmission probability and ABS 
peak height at the left and right interface. The subscript 
± corresponds to the sign of e^- The relative height of 
GAR and EG peaks in different magnetization geometries 
can be fully understood with this simple picture. 
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FIG. 5: (Color online) Sketch of the spin- up (positive axis) 
and spin-down (negative axis) Andreev states at the FM/SC 
interfaces for P (a) and AP (b) orientation of the contact 
magnetizations. CAR and EC processes are shown, and ex- 
pressions for obtaining their relative size are given in the table. 



bound states at isbi and ±£b2 (with \ehi\ > kb2|)- 
For the P case, similar expressions as in Fig. [5] hold, 
with Wp replaced by [^^;£W£]p, where £ G {1,2} labels 
the bound states. In the AP case, the EC contribution 
at zt^b^ is Tr[w£W£]Ap^ and the CAR contributions are 
asymmetric in voltage: for eVi, = {— ^bi, — ^b2, ^b2, ^bi} 
they are {-r'^ , -T'^j, -r^-/, -T'^}[wiWi]js^, with 7 = 
[w2W2/wiWi]^ [1^. The corresponding relative con- 
tributions to Gnl are {r, —Tj^rj, —T}{T — t)[wiWi]ap. 

In conclusion, we have studied subgap Andreev states 
in a FM/SC /EM setup. Due to an interaction between 
the bound states induced at the two SC/FM interfaces, 
their energetic positions depend on the relative interface 
magnetization orientation. This leads to marked asym- 
metries of the non-local conductance, which we explain 
in terms of CAR and EC processes with a simple picture 
based on the Andreev state positions and their weight, 
thereby clarifying the important role of Andreev states 
in non-local conductance experiments. 



As an example, we consider = I^^P = T and 

Kit/ = = ^ ^^^^ r < T < 1 in the P (AP) 

configuration sketched in Fig. [5l In the P case, there 
is a single spin-resolved state at energies ±£b, and ah 
states have the same weight Wp at the interfaces. The 
EC conductance peak thus scales as T^Wp (r^Wp) for 
positive (negative) energies, and the voltage asymmetry 
in the EC signal in Fig. HJ^a) is due to the different trans- 
missions for spin-up and spin-down. The CAR conduc- 
tance scales as — TrWp, both for eV^ = i^b, and is 
therefore symmetric in cVl. As a result, the total NL 
conductance in Fig. ID^c) changes sign so that one could 
switch between CAR or EC by tuning eV^. The split- 
ting of the conductance peaks for longer lengths [dashed 
curves in Figs. lU^a-c)] is due to the ABS repulsion dis- 
cussed in Fig. [3l^a). In the AP case, both a spin- up 
and spin-down state are present at every ABS energy 
so both CAR and EC contain the sum of two contri- 
butions: one from each spin-state. For EC the conduc- 
tance peak scales as 2T r [wW]ap^ while the CAR sig- 
nal scales as — (T^ +'?"^) ['^^]ap, where Wj^ (Wap) is the 
height of the smaller (larger) ABS peak (for Fig. IH solid 
lines, [^^;W]AP = 0.45Wp). These expressions are valid at 
cVl = i^b, so the NL conductance signals for i9r = — i^l 
in Figs.m^d-f) are all symmetric in voltage. In contrast to 
the P case, for the AP case the conductance peaks are re- 
duced when the distance between the contacts becomes 
larger (compare solid/dashed curves in Fig. |4]) because 
the weight wap decreases with the interface distance. 

In more general cases, where 7^ I'^^rI as in the dot- 
ted curves in Fig. HI the NL conductance contributions 
show peaks at four ABS energies (compare to Fig. [2]). 
Their relative weight can be determined analogously to 
the considerations above. E.g., for Fig. [2jc-d) there are 
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